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ABSTRACT 

We demonstrate that the current helicity observed in solar active regions traces the magnetic helicity 
of the large-scale dynamo generated held. We use an advanced 2D mean-held dynamo model with 
dynamo saturation based on the evolution of the magnetic helicity and algebraic quenching. For 
comparison, we also studied a more basic 2D mean-field dynamo model with simple algebraic alpha 
quenching only. Using these numerical models we obtained butterfly diagrams both for the small-scale 
current helicity and also for the large-scale magnetic helicity, and compared them with the butterfly 
diagram for the current helicity in active regions obtained from observations. This comparison shows 
that the current helicity of active regions, as estimated by —A • B evaluated at the depth from which 
the active region arises, resembles the observational data much better than the small-scale current 
helicity calculated directly from the helicity evolution equation. Here B and A are respectively the 
dynamo generated mean magnetic field and its vector potential. A theoretical interpretation of these 
results is given. 

Subject headings: magnetohydrodynamics (MHD) - Sun: magnetic helds - dynamo - interior - surface 
activity - sunspots 



1. INTRODUCTION 

The solar activity cycle is believed to be a manifes- 
tation of dynamo action which somewhere in the so- 
lar interior generates waves of quasi-stationary magnetic 
field propagating from middle latitude towards the so- 
lar equator ("dynamo waves" ) . The tradi tional expla- 
nation of this dynamo action (|Parkerlll955| ) is based on 
the joint action of differential rotation and mirror asym- 
metric convection which results in what has come to 
be known as the a-effect, base d on the helicity of the 
hydrodynami c convective flow (|Krause fc Radlerl Il980l: 
!MofFattlll978f ). This explanation is however not the only 
one discussed in the literature and, for example, merid- 
ional circulation is also suggested as an important co- 



factor of the a-effect, see e.g., IDikpati fc Gilmanl (l200l : 
IChoudhuri et~aT1 (pOOl . 

In turn, traditional dynamo scenarios based on dif- 
ferential rotation and the classical a-effect have to 
include a dynamo saturation mechanism. One of 
the most popular saturation mechanisms is based on 
a contr ibution to the a-eff ect from magnetic fluctu- 
ations ()Pouauet et all I1976D . A relevant quantifica- 
tion of this effect involv es consideration s of magneti c 
helicity evolution, e.g., Kleeorin etHI (|1995l l2003f ). 
A key role in the evolution of magnetic helicity is 
played by magnetic helicity fluxes dKleeorin et al.l [20001: 
Blackman fc Field] l2000at iBrandenburg fc Subramanianl 
2005T ). Of course, this scenario is not the only one 
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that ha s been suggested: for example IBlackman fc Field 
(2000bT): IBlackman fc Brandenburg! (|2003! ): IBrandenburg 
(2007) consider coronal-mass ejections as an impor- 
tant part of nonlin ear suppression of the dynamo, and 
iMitra et al.1 (1201 If) discuss the eff ects of losses via the 
solar wind. iChoudhuri et al.l (|2004l ) believe that the cur- 
rent helicity in solar active regions is substantially mod- 
ified when magnetic tubes rise up to the solar surface. 

A natural way to resolve such controversies is to de- 
termine relevant quantities such as the a-effect through 
observations, this providing a check on the various sce- 
narios. Such an option is now becoming realistic, start- 
ing from the 1990s when the first attempts to observe 
current helicity in solar active regions were undertaken 
(ISeehaferl 119901: IPevtsov et al.l[l99l iBao fc ZhandU99l 
Hagi no fc Sakuraill2004lL ~ 

Twenty years of continuous efforts by several obser- 
vational groups, with the main contribution coming 
from the Huairou Solar Station of China, has resulted 
(|Zhang et al.l [2010) in reconstruction of the current he- 
licity time-latitude (butterfly) diagram for one full solar 
magnetic cycle (1988-2005). From this butterfly diagram 
it is apparent that the current helicity is involved in the 
solar activity cycle and follows a polarity law comparable 
with the Hale polarity law for sunspots - but rather more 
complicated. In other words, dynamo generated mag- 
netic field is indeed mirror asymmetric and this mirror 
asymmetry is involved in the solar activity cycle and can 
be used to unders tand its nature (Klecorin eFaDHool 
iZhang et al.ll200l . 

What however needs clarification is which dynamo gov- 
erning parameter is traced by such a surface proxy as a 
measure of current helicity in solar active regions. A 
naive idea here is to identify this part of the surface 
current helicity with the current and magnetic helici- 
ties of the dynamo generated small-scale magnetic fields 
deep inside the Sun (say, in the solar overshoot layer), 
which suppress dynamo action. To start with some- 
thing definite, this n aive interpretation was applied by 
IKleeorin et ail (|2003h and lZhang et al.l (pOOl . 

Recent progress in observation which has resulted in 
butter lfy diagrams for t he current helicity in active re- 
gions (Zh ang et al.ll2010T) makes it possible to go further 
with this topic. The aim of this paper is to argue that 
the current helicity in solar active regions directly reflects 
the magnetic helicity of the large-scale dynamo generated 
field. We organize our arguments as follows. 

Based on the observed butterfly diagram for the cur- 
rent helicity, we here confront the naive interpretation 
with the available observations. We examine a mean- 
field dynamo model with dynamo suppression based on 
the magnetic helicity balance, obtain the correspond- 
ing current helicity butterfly diagrams, and also that 
of the large-scale magnetic helicity, and compare them 
with those observed. This comparison shows that the 
evolution of the large-scale magnetic helicity resembles 
the observational data much more closely than that of 
the current helicity of the small-scale fields. In order to 
demonstrate the robustness of this result, we also con- 
sider a more primitive dynamo model, with a simple al- 
gebraic a-quenching. From this model we calculate the 
large-scale magnetic helicity and compare it with the ob- 
servational butterfly diagram. We find that this fits ob- 
servations more or less as well as that from a dynamo 



model based on helicity conservation. We conclude that 
the major part of the observed current helicity in active 
regions is produced in the rise of magnetic loops to the 
solar surface. 

2. THE ROLE OF HELICITIES IN MAGNETIC FIELD 
EVOLUTION 

In the evolution of the magnetic field different helic- 
ities play different roles. Considering the small-scale 
velocity and magnetic fluctuations, u and b respec- 
tively, there are three helicities: (i) the kinetic helicity 
H u = (u-curlu) that determines the kinetic a effect; (ii) 
the current helicity H c — (b-curlb) that determines the 
magnetic part of the a effect; and (iii) the magnetic he- 
licity H m = (a-b), where b = curia. We use the angled 
brackets (...) to denote spatial integrals over all relevant 
turbulent fluid. These integrated helicities are used in 
the sense of average helicity densities. 

The total magnetic helicity, the sum Hm + H m of the 
magnetic helicities of the large- and small-scale fields, 
Hm = A-B and H m respectively, is conserved for very 
large magnetic Reynolds numbers. Here B — curl A is 
the large-scale magnetic field. Note that, on the con- 
trary, the current helicity H c is not conserved. On the 
other hand, the kinetic helicity H u is conserved only for 
very large fluid Reynolds numbers when the large-scale 
magnetic field B vanishes. The characteristic time for 
the decay of kinetic helicity is of the order of the turn- 
over time t = i/u of turbulent eddies in the energy- 
containing scale, £, of turbulence, while the character- 
istic time of the sm all magnetic helicity decay is of the 
order of T m = rRm dMoffattl[l9"78t iZeldovich et aH ll98l 
IBrandenburg: fc S^rb ramanianl l2005r ). where Rm = I u/rj 
is the magnetic Reynolds number, u is the characteristic 
turbulent velocity and rj is the magnetic diffusivity due 
to electrical conductivity of the fluid. The small-scale 
current helicity, H c is not an integral of motion and the 
characteristic time of H c varies from a short timescale, t, 
to much larger timescales. On the other hand, the char- 
acteristic decay times of the current helicity of large-scale 
field, He = BcurlB, and of the large-scale magnetic 
helicity, Hm, are of the order of the turbulent diffusion 
time. For weakly inhomogeneous turbulence the small- 
scale current helicity, H c , is proportional to the small- 
scale magnetic helicity, H m . 

As the dynamo amplifies the large-scale magnetic field, 
the large-scale magnetic helicity Hm = A-B grows in 
time (but not monotonically in a cyclic system). The 
evolution of the large-scale magnetic helicity Hm is de- 
termined by 

+ V • P M = 2£ ■ B - 2 V H C , (1) 

(IKleeorin et al.l [19951: IBlackman fc Field! I2000at 
IBrandenburg fc Subramanianl 12005ft . where £ = (uxb) 
is the mean electromotive force that determines gener- 
ation and dissipation of the large-scale magnetic field, 
25 ■ B is the source of the large-scale magnetic helicity 
due to the dynamo generated large-scale magnetic field, 
and Fm is the flux of large-scale magnetic helicity 
that determines its transport. Since the total magnetic 
helicity over all scales, Hm + H m integrated over the 
turbulent fluid, is conserved for very small magnetic dif- 
fusivity, the small-scale magnetic helicity changes during 
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the dynamo process, and its evolution is determined by 
the dynamic equation 



dH„ 
dt 



V F = -25 • B - 2rjH c , 



(2) 



(IKleeorin et all [l995t [Blackma n fc Field! [2000a; 
iBrandenburg fc Sub ramanian 2005), where —25 • B is 
the source of the small-scale magnetic helicity due to the 
dynamo generated large-scale magnetic field, F is the 
flux of small-scale magnetic helicity that determines its 
transport and 2r]H c = H m /T m is the dissipation rate of 
the small-scale magnetic helicity. It follows from Eq. (TTJ) 
and ((21) that the source of the small-scale and the large- 
scale magnetic helicities is located only in turbulent 
regions (i.e., in our case, in the solar convective zone). 
The magnetic part of the a effect is determined by the 
parameter \ c — T H c /{12np), and for weakly inhomo- 
geneous turbulence \ c is prop ortional to the magnetic 
helici t y: y c = H m / '(18-Krj^o) (IKleeorin fc Rogachcvskii 
119991: iBrandenburg fc Subramanian! 120051) . where p is 
the density and r] T is the turbulent magnetic diffusion. 

3. THE OBSERVED CURRENT HELICITY BUTTERFLY 
DIAGRAMS 

The observed butterfly diagrams of electric current he- 
licity for solar active regions during the l a st tw o solar 
cycles has been presented by iZhang et al.l ((2010) . The 
general structure can be described as follows. Current 
helicity is involved in the solar activity cycle and follows 
a polarity rule comparable to (however more complicated 
than) the polarity rule for toroidal magnetic field which 
in turn comes from the Hale polarity rule for sunspot 
groups. Migration of the helicity pattern is clearly visible 
and located near the toroidal field pattern. The wings 
of the helicity butterflies are slightly more inclined to 
the equator than the magnetic field wings, but the for- 
mer follow in general the latter. As a quantity quadratic 
in magnetic field, the current helicity in one hemisphere 
has the same sign for both cycles, with the opposite sign 
in the other hemisphere (a kind of unchanging dipolar 
symmetry). 

Fig. 1 shows the distribution of the average helical 
characteristics of the magnetic field in solar active re- 
gions in the form of a butterfly diagram (latitude-time) 
for 1988-2005 (which covers most of the 22nd and 23rd so- 
lar cycles). These results are inferred from photospheric 
vector magnetograms recorded at Huairou Solar Observ- 
ing Station. 

This, the longest available systematic dataset cover- 
ing the period of two solar cycles, comprises 6205 vector 
magnetograms of 984 solar active regions (most of the 
large solar active regions of both solar cycles). Of these, 
431 active regions belong to the 22nd solar cycle and 
553 to the 23rd. We have limited the latitudes of ac- 
tive regions to ±40° and most of them are below ±35°. 
The helicity values of the active regions have been av- 
eraged over latitude by intervals of 7° in solar latitude, 
and over overlapping two-year periods of time (i.e., 1988- 
1989, 1989-1990, 2004-2005). By this method of aver- 
aging we were able to group sets of at least 30 data points 
in order to make error estimations (computed as 95 % 
confidence intervals) reasonably small. In this sampling 
we find that 66% (63%) of active regions have negative 
(positive) mean current helicity in the northern (south- 



ern) hemisphere over the 22nd solar cycle and 58% (57%) 
in the 23rd solar cycle. 

There are some domains in the diagram where current 
helicity has the "wrong" sign with respect to the global 
polarity law. These domains of "wrong" sign are located 
at the very beginning and the very end of the wings. 

Con cerning alternative e xplanations of the observa- 
tions (IZhang et al.l I201C ). iMackav fc Van Ba llegooiicn 
(|2005|) : lYeates fc Mackavl (|2009ft describe current helic- 
ity in solar active regions in terms of generation of non- 
potential coronal structures by surface differential rota- 
tion. Note that the surface differential rotation cannot 
generate large-scale and small-scale magnetic helicity. It 
can only redistribute the existing magnetic helicity, as 
can any non-uniform large-scale motions. This process is 
determined by the flux term in the evolutionary equation 
for the magnetic helicity. The local change of the mag- 
netic helicity inside a given active region by the surface 
differential rotation plays an important role. In our pa- 
per we study robust global (rather than local) features of 
the evolution of the magnetic helicity, by averaging over 
an ensemble of active regions. In this case the global 
evolution of the magnetic helicity mainly depends on the 
sources of magnetic helicity inside the solar convective 
zone. 

4. AN ESTIMATE FOR THE CURRENT HELICITY IN 
ACTIVE REGIONS 

In this section we estimate the current helicity in 
active regions. There is a common belief that ac- 
tive regions are formed due to some nonaxisymmet- 
ric instability of ~ 100 kG magnetic fields in the 
tachocline, e.g., dGilman fc Dikpatl 120001: I C ally et all 
2003; P arfrevfc Menoul I2007D . However, the existence 
of such strong fiel ds and the role of t his mechanism re- 
main questionable (Br andenburgl2005l ) . Another promis- 
ing mechanism of formation of active regions is related to 
a negative effective magnetic pressure instability of the 
large-scale dynamo generated mag netic field. This insta- 
bility was predicted theoretica lly (jKleeorin et al.l [i~996l 
iRogachevskii fc Kleeorin I l2007|) and detected recently i n 



direct numerical simulations ( Brandenburg et "ail 120 111 ). 
The instability is caused by the suppression of turbulent 
hydromagnetic pressure by the mean magnetic field. At 
large Reynolds numbers and for sub-equipartition mag- 
netic fields, the resulting negative turbulent contribution 
can become so large that the effective mean magnetic 
pressure (the sum of turbulen t and non-turbulent con- 
tribu tions) appears negative ([Brandenburg et all I2010L 
1201 ID . In a stratified turbulent convection, this results 
in the excitation of a large-scale instability that results 
in the formation of large-scale inhomogeneous magnetic 
structures. This mechanism is consistent with the sug- 
gestion that acti ye regions are form ed near the surface of 
convective zone (|Kosovichevll2010[) . 

The spatial scale of an active region is much smaller 
than the solar radius, but much larger then the maximum 
scale of solar turbulent granulation. To estimate the cur- 
rent helicity in an active region, we have to relate the 
large-scale magnetic field B and its magnetic potential 
A inside the convective zone as well as the corresponding 
small-scale quantities inside the convective zone (which 
determine the small-scale magnetic fluctuations), with 
the surface magnetic field B ar and its magnetic poten- 
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Fig. 1. — Observed current helicity (white/black circles for positive/negative values) for solar active regions in the 22 nd and 23 rd solar 
cycles as averaged over two-year running windows over latitudinal bins of 7 degrees wide, overlaid with sunspot density (colour). The circle 
in the upper right corner of the panel indicates the typical value of observational uncertainty defined by 95% confidence intervals scaled to 
the same units as the circles. The vertical axis gives the latitude in degrees and the horizontal gives the time in years. 



tial A ar inside active regions, which are the quantities 
available to observations. 

We base our estimate on the following idea. Consider 
a new-born magnetic tube at the layer in convective zone 
from which the active region is developing. We assume 
that the rise of magnetic tubes is a fast adiabatic process. 
Let us also assume that the mean magnetic field and the 
total magnetic helicity vanish at the initial instant, and 
take into account the magnetic helicity conservation law 
(as solar plasma is highly conductive, and so we con- 
sider the magnetic helicity conservation law to hold at 
all scales, including that of the whole Sun). If this tube 
rises rapidly to the surface to produce an active region, 
the total magnetic helicity in the tube is conserved be- 
cause the process is rapid. Rising large-scale magnetic 
field and magnetic potential give the corresponding quan- 
tities for active regions, which may thus differ substan- 
tially from the corresponding quantities in the surround- 
ing medium. Because the initial total magnetic helicity 
of the tube, which was almost nonmagnetized, was neg- 
ligible, the magnetic helicity conservation law reads 

(A ar -B ar ) « -AB, (3) 

where the angular brackets denote averaging over the 
surface occupied by the active region. 

Now we relate the mean current helicity (B ar -curlB ar ) 
with the magnetic helicity (A ar -B ar ). We rewrite it from 
the first principles with the use of permutation tensors 
as 

(B ar -curlB ar ) ps -L (A ar -B ar ) + o(^f), (4) 

(see Appendix |A"]) . where R@ is the solar radius, and L ar 
is the spatial scale of an active region. Equations ([3} 
and (J2J) yield 

(B ar curlB ar ) a —L-A-B. (5) 



Therefore, the observed current helicity in active regions 
is expected to be a proxy for A • B. This idea will be 
checked using mean-field dynamo numerical modelling 
and comparison of the numerical results with the ob- 
served current helicity in active regions. 

5. DYNAMO MODELS 

Our approach to compare the dynamo models with ob- 
servations is as follows. We consider two types of dynamo 
models. Both types of models are 2D mean-field models 
with an axisymmetric magnetic field which depends on 
radius r and polar angle 9. The third (azimuthal) coor- 
dinate is (f>, and d/d<f> = 0. The dynamo action is based 
on differential rotation, with a rotation curve which re- 
sembles that of the solar convection zone, as known from 
helioseismological observations, and there is a conven- 
tional a-effect. 

The first type of model assumes a very naive algebraic 
a-quenching. Then we suppose that the total magnetic 
helicity is locally vanishing, so the magnetic helicity of 
the large-scale magnetic field produced in the course of 
mean-field dynamo action has to be compensated by 
small-scale magnetic helicity. (Thus we are assuming 
that at an initial instant the medium is non-magnetic, 
so that helicity conservation means that the sum of large 
and small-scale helicities remains zero.) We assume also 
that there is a separation of scales so that characteris- 
tic turbulence scales are much smaller than the charac- 
teristic spatial scales of mean magnetic field variations. 
This allows a link between current and magnetic helic- 
ities (jKleeorin fc Rogachevskiil [79991 to be made. This 
concept underlies the observational procedure for deter- 
mining the current helicity of active regions, and for cal- 
culating the current helicity from the magnetic helicity of 
the small-scale fields. Based on the same concept we es- 
timate the large-scale magnetic helicity as B^A^, where 

A(r, 6)4> is the magnetic potential for the poloidal field. 
As a result we obtain (for a given radius r) a theoretical 



Current and magnetic helicity of active regions 



5 



model for the current helicity as a function of t and 9 
which we overlay on the butterfly diagram for B^. We 
compare the result with the current helicity butterfly di- 
agram known from observations and obtained using sim- 
ilar underlying concepts. 

A further point is that the primitive model allows a 
simplification to the level of the one-dimensional Parker 
migratory dynamo, and t his opportunity h as been inves- 
tigated in this respect by Ku et al.l (|2009f ). We will use 
the results of that work for reference and comparison. 

We do not consider this primitive scheme as realistic. 
We are sure that any more or less realistic scenario for 
solar dynamo suppression will have to be much more 
sophisticated. On the other hand, we can see whether 
this primitive model produces a helicity butterfly dia- 
gram that is quite similar to that observed. The only 
shortcoming of the model is that the maximal current 
helicity occurs later then the maximum B<p, while it is 
observed to come up earlier. If magnetic helicity con- 
servation determines the nonlinear dynamo suppression, 
we expect that a careful reproduction of this balance, 
including helicity fluxes and the link between magnetic 
helicity and a-effect, will result in even a better theoret- 
ical butterfly diagram, and possibly improve the phase 
relations between helicity and toroidal magnetic field. 

As a specific example of the second type model that 
takes into account the influence of magnetic helicity 
balance on d ynamo action we use the dynamo model 
described by iZhang et al.1 (|2006| ) . Whereas simple a- 
quenching provides a quite robust suppression of a spher- 
ical dynamo and give (more or less) steady nonlinear 
magnetic field oscillations for a very wide range of pa- 
rameters, in contrast it is far from clear a priori that 
a dynamo suppression based on magnetic helicity con- 
servation is effective enough to suppress magnetic field 
growth and result in steady oscillations. In fact it works 
more-or-l ess satisfactorily only in a quite narrow p arame- 
ter range (jKleeorin et alj|2003t IZhang et al.ll2006l) . which 
appears inadequate to fit observations. 

We note two crucial points here. First of all, both types 
of models ignore any direct action of magnetic force on 
the rotation law. In the more primitive models, there is 
a crude parametrization of feedback onto the (purely hy- 
drodynamic) alpha effect [see Eq. ©] below. The latter 
more sophisticated model describes the back-reaction of 
the generated magnetic field on the dynamo process in 
terms of the magnetic contribution of the current helicity 
onto the magnetic part of a-effect. On the other hand, 
the feedback of the generated large-scale magnetic field 
on turbulent convection is described in our model by the 
algebraic quenching of a-effect, pumping velocities and 
turbulent magnetic diffusion. 

We assume that helicity conservation is not the only 
mechanism of dynamo suppression. The fact that we 
see a manifestation of helicity on the solar surface tells 
us that the buoyancy indeed plays some role, and we 
add it to the model. We stress that the buoyancy which 
we include in the model transports current helicity and 
magnetic helicity as well as large-scale magnetic field. 

Below we discuss the detailed dynamo models. We use 
spherical coordinates r, 9, (j) and describe an axisymmet- 
ric magnetic field by the azimuthal component of mag- 
netic field B 7 and the component A of the magnetic po- 
tential corresponding to the poloidal field. 



We measure length in units of the solar radius Rq , and 
time in units of a diffusion time based on the solar ra- 
dius and the reference turbulent magnetic diffusivity r] T0 . 
The magnetic field is measured in units of the equipar- 
tition field B cq — it* y/iirp*, the vector potential of the 
poloidal field A in units of RqB c(1 , the density p is nor- 
malized with respect to its value p* at the bottom of 
the convective zone, and the basic scales of the turbu- 
lent motions I and turbulent velocity u at the scale t are 
measured in units of their maximum values through the 
convective zone. The a-effect is measured in terms of 
ao, defined below, and angular velocity in units of the 
maximum surface value, f2o- 

5.1. The primitive, alpha- quenched model 
In the primitive dynamo model the a-effect is given by 



X v *v, 



(6) 



where \ v is proportional to the hydrodynamic helicity, 
H u , multiplied by the turbulent correlation time r, and 
= (1 + B 2 )^ 1 is the model for the a quenching non- 
linearity. For convenie nce, we use for most of these com- 
putations t h e cod e of iMoss fc Brook e (2000) - see also 
iMoss et ail (|2011l ). This code has the possibility of a 
modest reduction in the diffusivity, to ry m i n , in the in- 
nermost part of the computational shell (" tachocline" ) , 
below fractional radius 0.7. We define ?y r = i] ni i n /ri T0 . We 
also used this primitive formulation of alpha-quenching 
in the (otherwise very similar) model used in Sect. 15.21 
when producing Fig. In the latter case, the diffusivity 
is everywhere uniform. 

At the surface r = 1 the field is matched to a vacuum 
external field, and "overshoot" boundary conditions are 
used at the lower boundary. 

5.2. The model based on helicity balance 

Here we use the code described in IZhang et al. ( 2006), 
with two new features: we allow the possibility of merid- 
ional circulation and/or vertical motions attributed to 
magnetic buoyancy. The dynamo equations for A — 
r sin 9 A and B = r sin 9 B read 



dA 
dt ' 

dB 



{V 9 A + V e M )8A 



dt 



+ 



r 

d 2 A 
dr 2 
sm9 d 
r 



d6 
sin 9 d 
~d9 



89 



Vr~ = 

dr 
1 dA 
sin6> 39 
M )B 



C a aB 



(7) 



sin 9 



= C U 

d_ 

dr 



G 9 
&r d9 



r d 

or 



A 



d[V r B] 
dr 

sin 9 d 
' ~r 2 ~~d9 



Vb OB 



sin 9 d9 



8B 
dr 



(8) 



where V r = V r A + V r M + V B , V A (B) and V B (B) are the 
nonlinear drift velocities of poloidal and toroidal mean 
magnetic fields, V M is the meridional circulation ve- 
locity, Vb is the vertical buoyancy velocity, rj A (B) and 
r\ B (B) are the nonlinear turbulent magnetic diffusion co- 
efficients for the mean poloidal and toroidal magnetic 
fields, and and the non-dimensional dynamo parameters 
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are C a = a R Q /r] TO , C u 



The nonlinear 



turbulent magnetic diffusion coefficients and the nonlin- 
ear drift velocities are given in Appendix [B] The non- 
dimensional gradients of differential rotation are 



G r 



on 

dr 



on 
He' 



In this dynamo model with magnetic helicity evolution 
the total a-effect is given by 



a 



X 



(9) 



with a" — cto sin 8 cos 8 $„ . The magnetic part of the 
a-effect is based on the idea of magnetic helicity conser- 
vation and the link between current and magnetic helic- 
ities. Here x° an d X° are proportional to the hydrody- 
namic and current helicities multiplied by the turbulent 
correlation time, and $„ and $ m are quenching func- 
tions. The analytical form of the quenching functions 
& V (B) and $,„(£?) is given in Appendix IB1 The density 
profile is chosen in the form: 



p(z) — exp[— atan(0.457r,z)] 



(10) 



where z = I — — r) and /x = (1 — Rq/ Rq)^ 1 . Here 
a m 0.3 corresponds to a tenfold change of the density in 
the solar convective zone, a w 1 by a factor of about 10 3 . 



The equation for \ c = r sm & X° is 

2a 8A_8B_ 

2 ~m do 



Ox" 
8t 



T 



~sm8 8 ( 
'Va B — — I 



1 



1 

C a 
8A 



88 \sm8 88 

+iV A_ v B ) zdA + {v A 
or 



~8 2 A 
Ve K 88 



8 A 8B 
8r dr 



8{F r + (V B +V r M )x c } sin8 8 



dr 



08 



Fe + Vo M X c 
sin( 



(11) 

where T = r 2 sin 2 8T , T = F / (l8irr] T p) is related to the 
flux F of the small-scale magnetic helicity and given in 
Appendix [B] Rq /£ is the ratio of the solar radius to the 
basic scale of solar convection, T — (1/3) Rm (£/Rq) 2 is 
the dimensionless relaxation time of the magnetic helic- 
ity, Rm = £u/i] is the magnetic Reynolds number, with 
rj the 'basic' magnetic diffusion due to the electrical con- 
ductivity of the fluid. 

The meridional circulation (single cell in each hemi- 
sphere, poleward at surface) is determined by 



V, 



M 



1 



8[r^(r,8)} 



V, 



M . 



sin 8 r p(r) dr 
1 dV(r,6) 



sm8rp(r) 88 



(12) 
(13) 



R v sin 2 8 cos 8 f(r)p, f(r) = 2(r 
^2 



where ty(r,8) 

rb) 2 (r — 1)/(1 — r},) 2 , rj, is the base of the computa- 
tional shell. This is normalized so that the max of 
Vg 4 at the surface is unity. We introduce a coefficient 
R v — Rq Uo/r) TO , where Uq the maximum surface speed. 



Buoyancy is implemented by the introduction of a 
purely vertical velocity V B = in Eqs. (J7J, (JSJ 

and (HIl), where 7 > (|Moss et al. II1999H . We justify 
the apparent non-conservation of mass by adopting the 
argument of K.-H. Radler, pres ented as a private com- 
munication in lMoss et al.l (|1990) , that the return velocity 
will be in the form of a more-or-less uniform "rain". In 
some ways the process represents pumping by a " fountain 
flow" . As a result the regular velocity Vb appears in the 
governing equations for the large-scale magnetic field and 
magnetic and current helicities but not the the equation 
for density. From the viewpoint of probability theory, in 
the first case Vg is a mean quantity taken under the con- 
dition that an elementary volume is magnetized so it does 
not vanish, while in the second case this mean is taken 
without any condition and vanishes. In our opinion, this 
idea can also be constructive for other problems with 
magne tic helicity advective fluxes, e.g. iShukurov at al.l 
((200(1 . 

At the surface of the Sun, r = 1, we use vacuum bound- 
ary conditions on the field, i.e. B = and the poloidal 
field fits smoothly onto a potential external field. At 
the lower boundary (the bottom of the solar convective 
zone), r = ro = 0.64, B = B r = 0. At both r — ro and 
r = 1, we set 8x c /8r = 0, where \ c is proportional to 
the current helicity [see Eq. ([9])]. 

6. SIMULATED BUTTERFLY DIAGRAMS FOR CURRENT 
HELICITY 

We performed an extensive numerical investigation of 
the models in a parametric range which is considered to 
be adequate for solar dynamos. We estimate the val- 
ues of the governing parameters for different depths of 
the convec tive zone, using mode l s of th e solar c o nvect ive 
zone, e.g. iBaker fc Temesvarvl (fl966l ). iSpruitl (fl97l - 
more modern treatments make little difference to these 
estimates. In the upper part of the convective zone, say 
at depth (measured from the top) h* = 2 x 10 7 cm , 
the parameters are Rm = 10 5 , u = 9.4 x 10 4 cm s _1 , 
t = 2.6 x 10 7 cm, p — 4.5 x 10 -7 g cm~ 3 , the turbulent 
diffusivity r/ T = 0.8 x 10 12 cm 2 s _1 ; the equipartition 
mean magnetic field is B cq = 220 G and T — 5 x 10~ 3 . 
At depth /i* = 10 9 cm these values are Rm = 3 x 10 7 , 
u = 10 4 cm s -1 , t = 2.8 x 10 s cm, p = 5 x 10~ 4 g 
cm~ 3 , r\ T = 0.9 x 10 12 cm 2 s _1 ; the equipartition mean 
magnetic field is B cq = 800 G and T - 150. At the bot- 
tom of the convective zone, say at depth h* = 2 x 10 10 
cm, Rm = 2 x 10 9 , u = 2 x 'lO 3 cm s~\ I = 8 x 10 9 
cm, p — 2 x 10 _1 g cm~ 3 , 77^ = 5.3 x 10 12 cm 2 s _1 . 



Here the equipartition means magnetic field B c 



3000 



G and T w 10 7 . If we average the parameter T over 
the depth of t h e con vective zone, we obtain T w 5, see 
iKleeorin et all (pflol . 

We start with the results for the primitive model. 
Fig- El presents the current helicity butterfly diagrams 
overlaid on those for the toroidal field. We estimate this 
quantity based on the idea that the observed current he- 
licity in active regions is expected to trace —A • B, and 
so we plot in this section — A ■ B. 

We see that the plots successfully represent main fea- 
ture of the observed helicity patterns. The pattern pre- 
sented in Fig.[2]is quite typical for the model. Of course, 
one can choose a set of dynamo governing parameters 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Time, non dim units Time, non dim units 

Fig. 2. — Current helicity of active regions estimated as —A ■ B [see Eq. overlaid on toroidal field for the primitive dynamo model: 
left panel - deep layer, right - surface layer. C a = —6.5, C w = 6 ■ 10 4 , R v = (i.e. no meridional circulation), no buoyancy. The diffusivity 
constant ?j r = 0.5 and the bottom of the computational region is at ro = 0.64. The colour palette is hereafter chosen as follows: yellow is 
positive, red is negative, and green is zero. 



which is less similar to the observations. For example 
one can concentrate magnetic fields in the deep layer of 
convective zone (say, in the overshoot layer) by choosing 
a reduction rj r = 0.1 in the nominal "overshoot layer", 
instead of r] r — 0.5 (our standard case), as used in Fig.O 
This tends to make the helicity wave in the overshoot 
layer look more like a standing wave, but however keeps 
the main features of surface diagram (Fig.[3|). The highly 
anharmonic standing patterns of the butterfly diagram 
th at were discussed as a possible option for some stars, 
see lBaliunas et all (|2006f ) . look however irrelevant for the 
solar case. 

Of course, the helicity pattern in the butterfly dia- 
gram obtained in the models for particular choices of 
dynamo governing parameters can be slightly d i fferent 
from the observed helicity patterns. Ku et all (2009) 
demonstrated that meridional circulation can be used to 
make the simulated pattern resemble more closely that 
observed. 

We produced the same type of plots for models based 
on helicity conservation Fig. [5] We also present in Fig. @] 
the small-scale current helicity \c- Here the migration 
diagrams are presented for the middle radius of the dy- 
namo region: nearer the surface \ c displays only rela- 
tively weak vacillatory behaviour. We see that the small- 
scale current helicity is strongly concentrated in middle 
latitudes and helicity oscillations which are present in the 
model are almost invisible on the background of the in- 
tensive belt of constant helicity in middle latitudes. We 
doubt that such oscillations would be observable. We 
stress that, if this model produces any travelling helicity 
pattern, it is situated in the deep layers only. The pat- 
tern usually is much more similar to that presented in 
the right hand panel of Fig. [3] rather than to a travelling 
wave such as presented in Fig. [2] 

Taking our models as a whole, determination of further 
details of the location of the dynamo layer or the helicity 
sign inversion at some latitudes and phases of the solar 
cycle is beyond the scope of this investigation. Such a 
study would require a more complex simulation, which 
we hope to perform in the future. 

The numerical models used here are extremely simple 



and cannot be expected to sample the full range of so- 
lutions that are accessible to a three-dimensional, highly 
stratified, extremely turbulent system with complex and 
unknown boundary conditions. In the present paper we 
have only shown the theoretical possibility that the cur- 
rent helicity observed in solar active regions may trace 
the magnetic helicity of the large-scale dynamo gener- 
ated field. We demonstrated that the results obtained 
with our simplified dynamo model are compatible with 
the observations. At the same time we note that such 
computationally complex and expensive simulations are 
still far from being able to reproduce the range of ob- 
served solar phenomena, and that for the immediate fu- 
ture mean field models will continue to play an important 
role. Indeed, they may provide more immediate physical 
insight. 

The simple models which we have considered here were 
not intended to reproduce very fine details of the spatial- 
time distribution of helicity observable in the form of 
butterfly diagrams. We may again note that our mod- 
elling of helicity in the solar convective zone and active 
regions is still too simplified to be able to detect more 
detailed properties of the helicity dynamics. 

While we feel that it is remarkable that there is a class 
of kinematic models that are able to reproduce the cur- 
rent helicity observations from vector magnetograms, we 
may again note that some contribution to the observed 
current helicity may be due to the surface effects caused 
by the differential rotation (|Mackay fe Va n Ballcg ooiienl 
l2005HYeates fc Mackavll2009D and not from the dynamo. 

7. CONCLUSIONS 

We conclude that current helicity of magnetic field in 
active regions is a tracer of the magnetic helicity of the 
large-scale magnetic field in the solar interior. We be- 
lieve that this provides a unique option for tracing this 
quantity, which is very important fo r the solar dynamo . 
According to the observational data (Zh ang et al.ll2010l ). 
the current helicity in active regions is mainly negative in 
the Northern hemisphere. Numerical models give a nega- 
tive value for — A- B in the surface layer of the convective 
zone in the Northern hemisphere. 

Summarizing, we conclude that the current helicity 
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0.4 0.6 0.8 1 0.2 0.4 0.6 

Time, non dim units Time, non dim units 

Fig. 3. — Current helicity of active regions estimated as —A ■ B [see Eq. (0], overlaid on toroidal field for the primitive dynamo model 
with enhanced dynamo activity in the overshoot layer: left panel - deep layer, right - surface layer. The values of C a , C u , Rv and ro are 
the same as in the previous figure but diffusivity contrast r] r = 0.1. 




0.4 0.6 

Time, non dim units 



Fig. 4. — Small-scale current helicity Xc overlaid on toroidal field contours for the dynamo model based on helicity balance, near the 
middle radius of the dynamo region. C a = —5, Cui = 6 ■ 10 4 , R v = 10 (i.e. with meridional circulation), buoyancy parameter 7 = 1. The 
plots are for fractional radius 0.84. 
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of the magnetic field in active regions is expected to 
have the opposite sign to AB, evaluated at the depth 
at which the active region originates. Thus, the mod- 
els presented here are consistent with the interpretation 
that the mechanism responsible for the sign of the ob- 
served helic i ty op erates near the solar surface, cf. e.g., 
iKosovichevI (|2010l ). The mechanism of formation of the 
current helicity in active regions still requires further in- 
vestigation. 
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APPENDIX 

A. CURRENT HELICITY VERSUS THE MAGNETIC HELICITY 



Here we relate the mean current helicity (B al curl B ar ) with the magnetic helicity (A ar -B a 
mean current helicity from first principles with the use of permutation tensors as 



First we rewrite the 



(B ar curlB a 



z mm e mij lim V- Vf (A ar (x)5 ar (y)) = Um [(V x ■ V y ) (A ar (x) • B al '(y)) - V x V y (A«(x)B^(y))] : 



(Al) 

where £y„ is the fully antisymmetric Levi-Civita tensor, Rq is the solar radius, L ar is the spatial scale of an active 
region and r = x y. The use of the full tensor notation and limits are needed here in order to separate the large-scale 
and small-scale variables and to obtain a simple final answer in scalar form when there is separation of scales. Since 
r = x — y is a small-scale variable, R = (x + y)/2 is a large-scale variable, the derivatives 



d 

dXr, 



_d_ 



1 d 



2dR v 



-V y 

p 



d 



V y 

p 



d 

dy p 



1 d 



d 



dr v 2 dR v 



d 
dR v 



This implies that 
V x 



7V - 



dv 2 



i d 2 



4 9R 2 



V x V y = V y V x - V y — — 

v p v q v p v <3 v p Qft 



- V 



d 



d 2 



g dR n dR„dR„ 



We take into account that divB ar = [i.e. V y S ar (y) = 0] and divA ar = [i.e. V x A ar (x) = 0]. We also take into 
account that the characteristic scale of an active region is small compared with the thickness of the convection zone 
or the radius of the Sun. This implies that 



V x V y <^ r (x)S p ar (y)) 
and therefore this term vanishes. This yields 

(B ar -curlB ar ) = 



d 2 



8RpdR q 



(A ar (x)i? ar (y)) ~ O 



Rq 



dr p dr p 



(A ar -B a 



O 



4 

Rq 



(A2) 



Now we take into account that the second derivative of the correlation function 

d 2 

(A ar -B ar ) 



dr p dr p 



should be negative since as r — ¥ the correlation function has a maximum. Thus, we finally obtain 



(B ar curlB 3 



1 

L 2 r 



(A ar -B a 



O 



R 2 

-ft© 



(A3) 



Similar calculations relat i ng the current helicity and the magnetic helicity in k-space can be found in Appendix C of 
iKleeorin fc Rogachevskiil (| 1999( 1. 



B. DETAILED DESCRIPTION OF THE DYNAMO MODEL 
The flux of the small-scale magnetic helicity is chosen in the form 

T = 77, (B) B 2 {Cx V [ X v ^ (B)} + C 2 X v 4>v (B) A p } - C 3 k V X c , 
with A p = — Vp/p and k is the coefficient of turbulent diffusion of small-scale magnetic helicity (see below). 



(Bl) 
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The quenching functions <£>„ (B) and $ m (_B) appearing in the nonlinear a effect are given by 

3 
8B 2 



(B) + 3L(V8B)} 
arctan(v / 8i?) 



1 - 



V8B 



(B2) 
(B3) 



(|Rogachevskii fc Kleeorin~ll2000[) . where L(y) = 1 - 2y 2 + 2y 4 ln(l + y- 2 ). 

The nonlinear turbulent magnetic diffusion coefficients for the mean poloidal and toroidal magnetic fields, r) A (B) 
and rj B {B), are given in dimensionless form by 



r, A (B) = A 1 (4B)+A 2 (AB) , 

Vb (B) = A, (4fl) + ~ [2A 2 (45) - A 3 (45)] 



(B4) 
(B5) 



(|Rogachevskii fc Kleeorin"l l2004). For the case of weak magnetic field the turbulent diffusion coefficients are (in units 
of the reference value rj T0 ) 



Va =1-^B 2 , ^=l-325 2 , 



while for strong magnetic fields the scaling is 



Va = 



8B 2 



3\/2B 



(B6) 



(B7) 



The transition from one asymptotic form to the other can be thought of as occurring in the vicinity of B ~ J5 eq /4. 

The nonlinear drift velocities of poloidal and toroidal mean magnetic fields, V A (B) and V B (B), are given in 
dimensionless form by 



where 



V A (B) = V x (B) ^ + Wl (e r + C ot9e e ) + V p (B) 
2 r 



V 1 (B) = -A 3 (4B)-2A 2 (4B) , 
V 2 (B) = ^A 2 (4B) , 
V 3 (B) = ^[A 2 (4B)-A 3 (4B)} , 
V P (B) = Ia p [-5A 2 (4B) + Ia 3 (AB)} 



The asymptotic formulae for these velocities are given by 

As + 3A, 



V A = -B 2 



cot 9 eg 



for a weak magnetic field, and 



V 1 



V A = - 



32 



1 



B 



3A £ 



cot 9 eg 



V 1 



3VSB 

4 e. 



e r + cot 9 eg 



16B 2 p 



cot 9 eg 



(B8) 
(B9) 



3y/8B r l&B 2 p 

for strong fields. Here Ab = (VB 2 )/B 2 , e r and eg are unit vectors in the r and 9 directions of spherical polar 
coordinates, [A p ] r = —dlnp/dr, and [As] r = d\nB 2 /dr. See, for details, Rogachevskii & Kleeorin [2004, Eqs. (18), 
(19), (22)-(24)], which have been rewritten here in spherical geometry. 
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The functions Ak{y) are 



My)- 
MvY- 
Mv)-- 



arctan y 



6 

"5 
2 

V 2 



y 

arctan y 

y 

V 



15 



1 + vJ-7 L ^"v 



(y 2 + 3) - 3 



The nonlinear quenching of the turbulent magnetic diffusion of the magnetic helicity is given by 



k(B) 



l + A 1 {AB) + ^A 2 {AB) 



(BIO) 



The coefficient of turbulent diffusion of magnetic helicity K also has a dependence on B, namely k(B) = 1 — 24£? 2 /5 
for weak magnetic field and 



<B) = \\l 



3tt 
40B 



in the strong field limit. 
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